Current-oscillator correlation and Fano factor spectrum of quantum shuttle with finite bias voltage 

and temperature 



> 



O 



O 

o 



X 



Wenxi Lai, Yunshan Cao, and Zhongshui Ma 

School of Physics, Peking University, Beijing 100871, China 



fN| ' Abstract 

o 

(N 



A general master equation is derived to describe an electromechanical single-dot transistor in the Coulomb 
blockade regime. In the equation, Fermi distribution functions in the two leads are taken into account, which 
allows one to study the system as a function of bias voltage and temperature of the leads. Furthermore, we treat 



the coherent interaction mechanism between electron tunneling events and the dynamics of excited vibrational 
modes. Stationary solutions of the equation are numerically calculated. We show current through the oscillating 
island at low temperature appears step like characteristics as a function of the bias voltage and the steps depend 
on mean phonon number of the oscillator. At higher temperatures the current steps would disappear and this 
event is accompanied by the emergence of thermal noise of the charge transfer. When the system is mainly in 
the ground state, zero frequency Fano factor of current manifests sub-Poissonian noise and when the system 



CZ3 

is partially driven into its excited states it exhibits super-Poissonian noise. The difference in the current noise 
would almost be removed for the situation in which the dissipation rate of the oscillator is much larger than the 
bare tunneling rates of electrons. 



PACS numbers: 85.85.+j, 72.10.Bg, 73.23.Hk, 72.70.+m 



1. Introduction 



■ Nanomechanical oscillators has been the subject of active research during the past decade. They have potential applications 

■ in precision measurements [1] and quantum information |0. A typical example of the structures is electromechanical tunneling 
called quantum shuttle which transports electrons with the help of mechanically oscillating island 1 3 ] . Such devices with different 
sizes have been realized in experiments 10,0]. Theoretically, to derive a general equation for the electromechanical system is 

" appeared to be a difficult task, which in a way hinders to exploit some properties of the system. 

When the island is very small (within a few nm in diameter), the mechanical oscillator can be seen as a quantum system and 
I-V curves exhibit stepwise characteristics [4]. It is widely accepted in theory that the phenomenon can be interpreted based on 



the multi channels which are provided by the quantized modes of the oscillating island J^- llOH . In these theoretical studies the 
coherent coupling between electron transport and the vibrational modes, which might become important when the system works 
in the quantum regime, is not taken into account. 

In a fully quantum mechanical description of the electromechanical process ill ill , the coherent coupling is included with a cost 
of dealing with large matrix in the mathematical treatment. The quantum mechanical model of the coherent dynamics is further 



developed to investigate the shuttling mechanism [[12 



14H . According to the charge-position (momentum) correlation, motion of 
the quantum shuttle can be divided into the regimes of shuttling, tunneling and their coexistence fb?! . In the shuttling regime, 
electron transport is highly deterministic characterized by the extraordinary sub-Poissonian Fano factor [13]. Measuring the shot 
noise, transition between tunneling and shuttling can be identified J14I . As its extension to spintronics, based on the coupling 
between transport of spin-polarized electrons and mechanical degree of freedom of island the shuttle instability is predicted to 



appear two stationary domains such as vibronic and shuttling depending on the applied electric and magnetic fields 11511 . A good 
application of the quantum shuttle in spin-detection is predicted recently Jl^j . The above studies are carried out in the limit 
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1511 . or with the electron distribution functions which are independent of energy levels of the mechanical 



of large bias . 

oscillator Jl^lKSH . Further studies about the current fluctuation, especially at finite bias voltages and large range of temperatures, 
is still an important issue which are not adequately discussed in the past. 

In this paper, we will develop a very general master equation for the description of the electromechanical tunneling. We 
consider Fermi distribution functions of the electronic leads, in which discrete energy levels of the vibrational modes are in- 
volved. We will know later that when the bias voltage is not very large the distribution functions are sensitive to the oscillator 
levels. In addition, both the diagonal and off diagonal terms of the coupling between oscillator dynamics and electron transfer 
are incorporated in the equation and both of them are shown to be very important to describe the device. Considering position 
dependence of tunneling rates and all vibrational modes of the harmonic oscillator in the equation, we intend to overcome the 
shortcomings in a recent attempt of describing the electromechanical system where the position dependence is neglected and 
just two modes of the oscillator are taken into account [17]. Furthermore, we will show that our master equation can be obtained 
by treating the tunneling and the dissipation terms uniformly instead of with the approach of separated prescription adopted in 



the early derivations 111 



12 



17H . In the limit of low temperature and high bias voltage, the present master equation would be 
in accord with those given previously I12I, [14I1 . Using the equation, we found some important characteristics of the electrome- 
chanical tunneling. Current steps as a function of bias voltage is correlated with the averaged phonon number of the mechanical 
oscillator at low temperature. The steps are disappeared when temperature is increased to a certain quantity and this process 
is corresponding to the emergence of thermal noise. Voltage range of one step as a function of gate voltage is half of that as a 
function of bias voltage. At the low enough bias voltage, the system dominantly works in its ground state and zero frequency 
Fano factor appears to be sub-Poissonian. When the voltage is high enough to excite the system it manifests super-Poissonian 
noise. In the previously studied systems which restricted to the limit of large bias, the dot-lead couplings can be arbitrarily 
large. In the present work, the applied bias can be properly varied, but the practical studies are on the other hand limited to 
weak couplings with the low applied bias. In other words, the broadening of electronic level in the system due to the leads is not 
included. 

Left lead , . , Right lead 




FIG. 1: The model illustration. A single-level dot is in the harmonic potential between two leads with bias voltage V . One possible path of an 
electron is marked with the arrows. The dashed lines imply trajectory of the electron, and the dotted line denotes damping of the vibrational 
mode. 



2. The system model and derivation of the master equation 

We consider a model that a single-level quantum dot (QD) is connected to two leads by some elastic molecules. Mechanical 
vibration of the dot is described by a harmonic oscillator with an effective mass m. With bias voltage, V, between the two 
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leads, electrons can be transferred from one lead into another. The sketch of the model is shown in figure Q] With respect to the 
energy level of the QD, the chemical potentials in the left and right leads are eV/2 and —eV/2, respectively, where e is absolute 
value of the electron charge. We assume the capacitance of the QD is so small that its electron occupation number is or 1. 
Initially, the QD is empty and the harmonic oscillator is in its ground state. When an electron jumps onto the QD, the electric 
field caused by the bias voltage exerts a force on the charged dot and drives the mechanical oscillator into its excited states. As 
a consequence, the electron transfer is influenced by the vibration. The total Hamiltonian of the model is sum of the electron 
tunneling Hamiltonian Htun, the mechanical oscillator Hamiltonian H mec h and the coupling between charge and the oscillator 

Hdriv - 

H = H tun + Hmech + Hdriv (1) 

The electron tunneling Hamiltonian is written in the form 

Htun = £oC f e + J2 Zykdlkdyk + fr ( T yke S » a{af+a) dl k c + T; k e s y a ^+^cU yk ), (2) 

k\y — l,r k;y—l,r 

where the first term is energy of the QD with annihilation and creation operators c, cK The second term describes noninteracting 
electrons in the left (y — I) and right (y=r) leads. The operator d yk {d yk ) creates (annihilates) an electron with momentum fc in 
the lead y. Spin degree of freedom is not involved here. The third term represents electron tunneling between the leads and the 
QD. The coupling strength is exponentially depend on the dot position and the operator of coordinate is denoted in the second 
quantization form with a and aJ. For simplicity, we define Si = —l,S r = 1 and a = ^ which is fraction of the zero point 
position uncertainty xq = 2w ^ o and the tunneling length A. The mechanical oscillator Hamiltonian is described by 

H m ech = hw a^a + 2J frio k b\b k + 2J hg k {b\a + a)b k ), (3) 

k k 

where the first term is free evolution of the mechanical oscillator with the inherent frequency oj . The creation and annihilation 
operators b\., b k in the second term characterize the Bosonic thermal bath. The last term represents oscillator dissipation due to 
the bath. The localized charge in the QD is coupled to the oscillator in the form 

H dHv = -ftfi(a t + a)c t c, (4) 

where fl = eV h *° and d is effective distance between the two leads. For simplicity, we assume only bias voltage V contribute to 
the electric field. Other electric environment of the QD may cause bias voltage independent electric field. It also contributes to 
the charge-oscillator coupling, which mainly plays the role of gate voltage by shifting the QD level and changing the number of 
vibrational states accessible for the electron transport We are not interested in this effect and so discard this term here. 
Now, to derive the master equation, we rewrite the Hamiltonian into a free Ho and an interaction part Hdriv + H\, where 

H =e c^c+ ^ jyhd] lk dyk + fatWa + tjLJ k blb k , (5) 

k;y—l,r k 

and 

H i = H — H — Hdriv ■ (6) 

In the interaction picture, the total density matrix pr(t) is pr(t) = e ltH °/ h pT{t)e ~ ltH o/ h which satisfies the Liouville- 
vonNeumann equation 

= ±.[H driv (t),p T (t)} + (7) 
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where H drlv (t) = e aH °/ h H drlv (t)e- UHo/h and H x {t) = e ltHo/h H 1 (t)e- ltH °/ h . Integrating equation (7) overtime t-y and 
substituting it into the second commutator of equation (7), we obtain 

= ^[H driv (t),Mt)} + ^[Hi(t),p T (0)}+(^) f [ffi(t),[ff d ri,(ti)+-ffi(ti),Pr(ti)]]<iti (8) 
at in in V"/ Jo 

It is pointed out specially that equation (8) is deduced in a way which is a bit different from the the standard approach for the 
derivation of master equation [23]. The first commutator in the right side of equation (7) is retained in equation (8) considering 
linear response of the motion of charged dot to the electric field. Meanwhile, the equation involves higher order of Hi to describe 
interaction between the opened system and the large reservoirs with a great number of degree of freedom. 

There are two kinds of environment for the system, one is electronic reservoir as the two leads of electrons, another is the 
Bosonic bath which is interacting with the mechanical oscillator. Under the Born approximation, the total density matrix can be 
factorized as pr(t) = p(t)pLPB, where p(t) is density matrix of the system composed of the QD and the mechanical oscillator. 
Both of the leads and the Bosonic bath are assumed to be in equilibrium state all the time and descried by the density matrix pl 
and pb, respectively. Therefore, if initially pt(0) = p(0)plPb, then at time t, pr{t) = p{t)PLPB- This approximation is valid 
as long as the coupling is weak and the environment is so large that the back action from the system to the environment can be 
negligible. Performing trace over the leads (tr^) and bath (trs) variables we obtain the reduced density matrix for the system 
p(t) = trBtrL{pT{t)}. Besides, the reservoirs are assumed to have much short correlation time compared to the coherent 
evolution of the system. Hence, there is no memory about the system history available in the reservoirs. As a result, the system 
evolution does not have to do with its history. Then, in the Markov approximation, we replace p(ti) by p(t), and arrive at the 
equation for the reduced density matrix 



= ^[H d riv{t),p{t)]- tTBtTL[Hl{t) t [Hi{ti),p{t)p L p B \\dtx (9) 

After evaluating the commutation relation of equation (9), we transform it back to Schrodinger picture using the unitary 
operator e -**( E oc t c+fiw ci t o)/ft i p or convenience, we make a time displacement as r = t — t±. Denoting electron number in 
the left and right leads as w and v, respectively, we write the density matrix in the form p w,v (t). If an electron is annihilated 
(created) in the left lead it becomes p w ~ 1 < v (t) (p w+1,v (t)), and if an electron is created (annihilated) in the right lead it would 
be p w ' v+1 (t) (p w ' v ~ 1 (t)). Additionally, the Fermi distribution function in lead y is introduced by tr^ld^dykPL} — fy{£yk), 
where fy(£ y k) — l/( ex P[/3(?j/fc — Py)] + 1) with y = l,r and /3 — 1/fcsT. The Bose distribution function in the thermal bath 
of mechanical oscillator is obtained from trs{b\bkPB} = n>B(^k), where nsi^k) = l/(exp[/3fio;fe] — 1). Then we reach the 
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following equation 



f)n w ' v (t) 1 
P K ' = —[e Jc+hwo^a~hn(a^+a)Jc,p w ' v (t)] 



+ f dr V \T yk \ 2 e a ' 

JO u.,.-i „ 



/2 



(S y a) 



mi +ni 



k\y—l,r mi,ni— 



x [/ y (^ fe )e s « Q ( at+a )c t p , " + ^ 1_5y ^ 2 ' 1 ' + ^ 1+s ' ! '^ / ' 2 (t)c(a' t )" 11 ( a )"i e ~ l (5»fe^ e o+( m i^"i)' iLl 'ti) r /' i 

+ (1 - / y (^j /fc ))e 5yQ ^ at+a ^cp u,_ ^ l_s ' ! '^ 2 ' , ' _ ^ :L+S! '^ 2 (t)c 1 '(a' l ') mi ( a ) n i e i (S«'=- e "-( m i^ ,l i)''^o)T/?i 

+ (1 - fy{£,yk))(^) mi (a)" 1 cp tt ' _ ^ 1_5y '/ 2 '' u_(1+5y ^ 2 (t)c t e 5yQ ^ at+a ^e~ i ^ yfc_eo+(mi_ni ^" ' |r / ;i 

-Afefe)e 5yQ(at+a) (a t ) mi (a)™ 1 cc t p u '' l 'We l(5! " I ^ ~ (mi ^" l)n " o)T/n 

-Afefe)p 1 "' 1, Wcc t (a t ) mi (a) ni e SyQ(at+a) e" i(iyfc " eo+(mi " ni) ^ o)T/;i 

-(! - /y(^fe))e SyQ(at+Q Ha t ) mi (a)™ 1 c t cp tt, ' t '(i)e~ l( « yfc ^ 0+(mi -™ l) ^ o)T/?i 

-(1 - /. y (^ fe ))p w ' u (i)c t c(a t ) mi (a)™ 1 e 5yQ(at+a) e i( « yfc ~ eo ~ (mi ~" l)?i " o)r/R ] 

Jo k 

+ (1 + n B {uj k )){ap w - v {t)a)e l ^ k -^ o)T - a) ap w > v {t)e- l{Uk -" a)T ) + H.c] (10) 



Considering rapid decaying electrons in the reservoirs, the time integral regime in the equation is extended to infinite, i.e. 
t — > oo. The integrals over time r are performed using the formula J °° d,Te ±%XT = ±*'-Pj + 7r<5(x) for any variable The 
imaginary part ±iP- is corresponding to Lamb shift in quantum optics [23]. Since we consider weak coupling between the 
system and its reservoirs (electronic leads and Bosonic thermal bath), the shift is very small and neglected. Subsequently, by 
defining the density of states in lead y as N y (£ y k) and density of states in the thermal bath as D(ujk), we make the conversion, 
12k-y=i r ~ / d^ykNy^yk), k ~ I k D (ui k) , and integrate over the variables of electron and phonon, respectively. Under 
the wide band approximation, the tunneling rates is written as T y — 2tt \T y \ 2 N y (y=l,r) and the dissipation rate is 7 = 2nDg 2 , 
which are independent of energy spectrum. Finally, we achieve the master equation: 

dp 1 

— — — — [eqc'c + hu!Qa)a — ft£l{a) + a)c^c, p v ] 
at in 

1 „ . . (q n \m 1 +n 1 +m 2 + n 2 

+ - e a r V [ y ] 

2 v TOi!ni!m 2 !n 2 ! 

y—l,r mi,ni,m2,n2— 
x\f <A + n v + ( 1 + S v)/ 2 A- _ A- A+ n v 

+A+ p v+^+s y )/2 A - -p v A- A+ ) 

'l 1 ■) y,mim J\- e± mimf J ■™-m 2 n 2 - tl m 2 n 2 m 1 n 1 l J 

+A~ p»-{i+s v )/ 2 A+ -p»A + A- )] 

1 m^n^r JA nimi r A *~ri\m\ mini J J 

+n Bl D[a)]p v + (1 + n B ftD[a]p v (11) 

where, without loss of generality, we just denote the number of electrons v collected in the right lead, which is required in the 
present paper. The commutator on the right side of the equation denotes the system evolution which is driven by an electric 
field. The second part represents phonon assisted tunneling, in which the operators A^ n and n are defined as j4~ n = 

c(a') m * (a) riz , = c^(a^) mz (a) n *>, where z = 1, 2. A^ n means an electron is out from the dot along with m z phonon 
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are created and n z phonon are annihilated, and so on , The Fermi distribution function is written in the form fy >miTll — 

e $le +(m 1 - ni jhu, +Syev/2\ +1 ■ The last two terms are damping of the harmonic oscillator due to its coupling to the thermal bath 
with the Bose distribution function n B = l/(e phu ° - 1). D[a]p v is defined as D[a]p v = ap v {t)a) - \{a*ap v {t) + p" (t)a*a). In 
the present work, we take the same temperature for the electronic leads and the thermal bath. It is easy to extend to the situation 
in which temperature in the two kinds of reservoirs are not the same. 

3. Quantized current 

In this section we study the stationary characteristics of electron transport across the shuttle junction. In detail, the correlation 
between current and energy of the vibrational modes is exploited. Current with respect to the gate voltage is described which is 
good comparison with the I-V curves. We will analyze the consequence that arising from the dissipation of mechanical oscillator 
and temperature of environment which is composed of electronic leads and thermal bath of the oscillator. 

Because of the charge conservation, the stationary current can be calculated from the flow either in the left or the right lead. 
For convenience, we consider electrons in and out from the right lead. The probability of v electrons collected in the right 
lead is P" = tr mec h[tr c h a r[p v ]], where tr mec h is trace over the variables of the mechanical oscillator, while tr c h ar denotes 
trace over the charge degree of freedom in the QD. According to the counting theory Ivh . the current is achieved using the 
formula / = e Y^vLo V P" ■ Due to large Hilbert space of the fully quantized system considered here, analytical expression for 
the current vs voltage characteristics is impossible to derive, even in an approximate form. What we deduce from the current 
formula is an expression in terms of the density matrix. The density matrix elements are obtained by solving the stationary 
solutions of equation (11). To solve the equation, we project it onto the Fock state bases spanning the Hilbert space of the 



system ||24j]. The state bases of the Hilbert space are {\0} char , |l} c hor} 0{|O) mec h, \l}mech, —\ n )-mech, •■•}• Where, (g) means 
direct product, 1 0} c har and 1 1 ) c h ar are eigenstates of the QD, indicating the dot is occupied by zero and one electron, respectively, 
and \n) mec h is eigenstate of the ?ith level of the mechanical oscillator. In this state space, the density matrix elements are denoted 
as poo, mn = m ech (m\ char (0 1 p\ 0) char | n) mech for empty dot and pu >mn = mec h {m\ char (l\p\l) char \n) mech for the dot occupied 
by one electron, where m, n — 0, 1,2, ... Obviously, the density matrix implies the coupling between electron hoping and 
coherent motion of center of mass. After the projection, assuming dp/dt = 0, we would reach a system of linear equations of 
the density matrix elements. Associated with the normalization condition X)^Lo ^ = 1> we numerically solve the system of 
linear equations directly. The effective dimension of the Hilbert space is directly related to the excited states of the harmonic 
oscillator which are involved in practical transport. In the present paper, we consider 15 lowest vibrational modes for the 
numerical calculation, that is < to, n < 15. It is enough for the system which works under the low bias voltage and weak 
dot-lead couplings. For the parameters adopted here, the probability spectrum of phonon distribution, poo >nn + Pn.nn, implies 
that contribution from the higher levels (n > 15) of the vibrational modes is negligible small. Actually, the master equation 
given in this paper is valid for much wide parameter range than that used here as long as adequate number of vibrational states 
are taken into account. However, the more number of vibrational modes are considered, the longer time is cost in the numerical 
implementation. The large memory and long time requirements are weak points of the approach to directly solve the equation. 
These problems are believed to be circumvented in the iteration method in which a preconditioning is necessary to ensure the 
convergence iHl . The iteration reaches its end when sum of the diagonal elements of the system density matrix closes to unity. 

The current as a function of bias voltage is plotted in figure [2] (a) with the solid line. Almost discontinuous transitions in 
the low temperature are found. It is shown that no current is available at zero bias voltage. When the bias moves away from 
the zero point to a small quantity, current appears and sharply increases. And then, even the voltage keeps rising, the current 
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FIG. 2: (a) J — V curves of the system is illustrated. The solid line is current and the dotted line is corresponding differential conductance, 
(b) The solid line is average number of phonon vs bias voltage. The dotted line is absorbtion (or emission) peaks of energy of the harmonic 
oscillator. Parameters in the two figures are the same, and they are T; = T r — O.OOlwo, xg/d = 0.003, a = 0.75, 7 = 0.02a;o, eo — and 
fifrwo = 20. 

seams to be unchanged. The chemical potential in the left lead now lies between the ground and the first excited states of the 
mechanical oscillator (see figureQ]). When the bias voltage rise to 2tkOo, the chemical potential reaches the first excited level of 
the oscillator, which results second jump of the current. More current steps emerge for further increase of the bias voltage. Each 
time the chemical potential reaches an additional level an extra transport channel is opened and a current step would be observed. 
The energy spacing of the oscillator levels is reflected by the voltage range of the current steps. Hight of the steps tend to be 
more and more small with the increase of the bias voltage and becomes invisible. In fact the hight of the steps can be controlled 
by tuning the parameters Ti, T r , 7 and a. The dotted line represents differential conductance corresponding to the current. The 
multi vibrational modes result in many resonant conductance peaks in quantum shuttle. Every peak implies a sharp increase 
of the current, between two peaks the current is stationary with respect to the bias. As the left and right parts of the model is 
absolutely symmetry, the current is appeared to be antisymmetry for the positive and the negative bias. The antisymmetry would 
be broken if the left (T;) and right (T r ) bare tunneling rates are not equal 

To further understand mechanism of the quantum shuttle, we move our attention to the mechanical oscillator. The key of 
the theoretical model is that the mechanical oscillator introduces multi modes into the dot conductor and quantizes the current. 
Discrete levels of the oscillator play the role of multi sub-bands in narrow conductor of two-dimensional electron gas in which 
the quantized conductance are observed [18]. Figure|2](b) shows averaged phonon number of the harmonic oscillator, which can 
be calculated from the formula (n) = X^^Lo n (Poo,nn + Pn,nn)- This curve corresponds to the stepped current in figure|2](a) 
with the same parameters. It is clear that every emission (absorption) of energy by the oscillator is accompanied by an increase 
(decrease) of current in the system. Therefore, the current is intensively correlated with energy of the mechanical oscillator. 
Similar relationship of the correspondence is discussed quite recently for asymmetric coupling between the island and the two 
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electronic leads 12011 . In the voltage area of — 2tiwo < eV < 2Hu!q, the mechanical oscillator is mainly in the ground state 
with zero mean phonon number, but it still contributes to the current due to the zero point fluctuation. We can easily prove it in 
figure |2] (a) that the first step is larger than the current of bare tunneling TiT r / (T; + T r ). The resonant absorption and emission 



peaks of the oscillator are illustrated by the dotted line in the figure. 




FIG. 3: The current vs energy level of the QD at bias voltages eV = 5fiuo (solid line), eV = 7?kJo (dashed line), and eV = 9hujo (dotted 
line). The other parameters are Ti = r r = 0.005wq, Xo/d = 0.003, a = 0.75, 7 = 0.03cjo and fthuo = 20. 

In the above analysis, we set the energy level of the QD as a constant. Now, we fix bias voltage and vary the level of the QD. 
In experiment, it can be realized by tuning the gate voltage. We plot the current as a function of the dot level in figure [3] for 
different bias voltages. As shown in the figure, area of the transport window is equal to the bias voltage eV . Outside the window, 
the current reduces to zero and the electron transfer is prohibited. Inside the window, we obtain current with stepwise structure. 
The highest step corresponds to the QD level that located in the middle of the bias. One need to shift the gate voltage by Hujo to 
obtain a new current step instead of 2Hljo in the bias voltage (see figure^. The figure implies that, by adjusting the gate voltage, 
one can control the current with discrete quantities. 

Current of the system is correlated with energy of the mechanical oscillator as discussed in the preceding subsection. Damping 
of the vibrational mode influences mean energy of the oscillator and so the electron transfer. Figured] illustrates such effect in 
the current by changing the dissipation rate. At the bias voltages eV = Stkvo, 5tkOo, and 7ftioo, the system is driven by the 
electric field and excited into high levels. Current across the excited system is remarkably decreased when the dissipation rate 
rises from a small quantity to about 10T; (and 10T r ). When the rate is higher than around 20T/ (and 20T r ), the currents 
approach a stationary value. At eV = Hloq, the mechanical oscillator is mainly in the ground state, therefore, the current is 
almost independent to dissipation rate (see the solid line). Approximate calculation reveals the currents for 7 < O.Olwo are 
further intensified but have finite quantities. 

Next, we consider temperature dependence of the system. As illustrated in figure|5](a), critical low temperature is required to 
observe the quantized current. If we increase temperature the step like structure tends to be smeared. Especially, the plateaus 
are disappeared at the temperature (3hu>o — 3. In the regime of high bias in this figure, the smoothed line (the dot-dashed line) 
is above the other curves that at the lower temperatures. In fact, it is not always the case. With high bare tunneling rate (e.g. 
Ti = T r = O.Olwo), the smoothed current would be lower than the stepped ones (not shown). Such temperature dependent 
effect is previously studied both experimentally I21I and theoretically |22]. Nevertheless, it makes sense to show it as we will 
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FIG. 4: (Color on line) The current vs the dissipation rate. Current across the excited system is denoted by the dotted, dashed and dot-dashed 
lines. Current across the system in the ground state is represented by the solid line (T; = T r = 0.005o;o, xo/d = 0.003, a — 0.75, eo = 0, 
/3hu = 20). 




FIG. 5: (a) (Color on line) I-V curves at different temperatures. With the increase of temperature, the stairs are blurred and disappeared in 
the end. (b) Zero frequency Fano factor as a function of temperature. The other parameters in the two figures are the same, i.e. T; = r r = 
O.OOl^o, x /d = 0.003, a = 0.75, e = and 7 = 0.02o;o- 

connect it to the behavior of zero frequency current noise in the next section. The step disappearance can also be seen in the case 
where the frequency-independent quality factor of the vibrational mode is very small 

4. Current fluctuation 
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FIG. 6: Fano factor spectrum at the dissipation rate 7 = 0.001cj - (a) F(0) < 1, (b) F(0) < 1, (c) F(0) > 1, and (d) F(0) > 1. The rest 
parameters are Ti = T r = 0.005o;o, :r /d = 0.003, a = 0.75, e = and ,9?ia;o = 20. 
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FIG. 7: Fano factor spectrum at the dissipation rate 7 = 0.3oj . (a) F(0) < 1, (b) F(0) < 1, (c) F(0) < 1, and (d) F(0) < 1. The rest 
parameters are the same as that in Fig. [6] 



In this section, let us focus our attention on the Fano factor, F(ut) = S(uj)/2el, which is given by ratio of the actual noise 
spectrum if? (w) and the Poissonian noise 2el. Concerning to the probability of v electrons collected in the right lead, we employ 
the McDonald formula UM, S(w) = 2e 2 w / °° dtsin(ujt) ^ Y,7=o v 2 P v {t), to calculate shot noise of the system. Figure [6] 
shows Fano factor spectrum for the dissipation rate which satisfies 7 < T/, T r . At the bias voltage eV — Hujq, mainly the 
channel of ground state is opened and the probability of electrons passing through the system is concentrated in this channel. 
Hence, the zero frequency Fano factor is appeared to be sub-Poissonian as illustrated in figure[6](a). It is similar to bare tunneling 



11 



process. The difference is that noise peaks appear at frequencies ±uq. It indicates an additional channel with small probability 
contributes to the transfer. Figure|6](b) reveals that, at voltage eV — 2fru)o, Fano factor is appeared to be super-Poissonian at the 
frequencies which are integral of cjq. At the voltage, excited states of the mechanical oscillator begin to effectively contribute 
to the electron flow, which destroys the zero frequency sub-Poissonian statistics. This effect can be observed more clearly at 
higher voltages(see figure [6] (c) and (d)). Since in the higher voltages, more channels are opened for the transport. The super- 
Poissonian noise at the zero frequency implies that the electron transfer through the ground state channel is interrupted by the 
tunneling through the channels of the excited levels. Due to the Coulomb blockade effect, there will be a competition between 
these channels with different tunneling probabilities and correlation is occurred among these transport processes of the multi 
channels. Electron transport in the present situation is not so deterministic as that in the so called shuttling regime where zero 
frequency sub-Poissonian noise is predicted It may be due to the extremely large bias voltage applied in their model, 

consequently, multi channels are not involved and electrons are forced to transfer mainly in one direction without reflected back 
to the original lead. In the area of off resonant frequencies, we can see noise suppression which is consistent with the result 
achieved from the model of incoherent dynamics B27I1 . 

In figure [7] Fano factor is shown for the situation 7 3> Ti, T r . Due to the fast damping effect, the contribution from excited 
states of the mechanical resonator is very small. A electron transports with the dominant probability through the channel provided 
by the ground state of the system. It causes suppression of the noise. Especially, the zero frequency Fano factor exhibits sub- 
Poissonian approaching 0.5 even in the case that the bias voltage is increased to a finite quantity. In fact, the result is also held 
in the limit of large bias voltage fl^ll . 

The above interpretation on physical picture of the super-Poissonian statistics is consistent with the previous results. In 
a movable QD array I25I . the different current channels are formed due to different resonant quantum states connecting the 
neighboring dots in the co-tunneling regime. The switching between those channels gives rise to super-Poissonian noise in the 
regime of small damping rate. In the semiclassical case, electron transport through a bistable coexistent channels of shutting 



and tunneling causes super-Poissonian noise spectrum both at zero [13] and finite frequencies [28]. The bistable of the quantum 



shuttle is further illustrated with Full counting statistics M29I1 



Now, let us briefly discuss the influence of temperature on the zero frequency current fluctuation. As illustrated in figure[5](b), 
large Fano factor is predicted in the case fihujQ ~ 1. The curves at different bias voltages have a common feature that they do not 
obviously depend on temperature until (3fkuo is decreased to about 3. When (3fkuo < 3, the noise is supposed to be dominated 
by thermal noise. We connect the Fano factor with the I-V curves at different temperatures in figure [5] (a). It is revealed that 
noise increase with temperature is accompanied by disappearance of the current steps near /3tki>o = 3. It seams to imply that the 
thermal noise which emerges due to the finite temperature removes the quantum mechanical characteristics of the current. 

5. Discussion and conclusions 

First we discuss the effect that comes from the coherent coupling between the charge transport and dynamics of the mechanical 
resonator. Now we calculate the system current using an incoherent model in which only diagonal terms of the density matrix, 
such as poo,nn, Pxx,nn, m the equation (11) and the current formula are involved. As illustrated in figure [8] for the lower bias 
eV < ifaujo, contribution from the off diagonal part is negligible small. However, in the case eV > 4:hajQ, the current calculated 
from the coherent model including both diagonal and off diagonal terms of the density matrix is lower than that achieved from the 
incoherent model. The suppression of current in the coherent model may be due to the destructive interference between different 
transport channels. The incoherent model applied here is not absolutely the same as those considered previously because of the 
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FIG. 8: (Color on line) Comparison between the currents calculated from the coherent and incoherent models. The solid and dotted lines 
present the currents which only involve diagonal elements of the system density matrix. The dashed and dot-dashed lines indicate the currents 
that calculated from both the diagonal and off diagonal terms of the density matrix. Here we denote T; = r r = T. The rest parameters are 
x /d = 0.003, a = 0.75, e = 0, and /3/kj = 30. 



different derivation methods |6- 10]. Whereas, here we just intend to make clear the importance of the coherent coupling in the 
system. 

The relation between resonator state and charge transport is an interesting issue. Measuring the character of charge transport 
in the system one expects to find out the information about mechanical resonator. In sections 3, we show that the current is 
sensitive to the mean phonon number of the resonator with the varying bias voltage. Then, we point out that the Fano factor 
spectrum of the charge transfer is dependent on the motion of mechanical oscillator. Recently, the current noise together with 
the phonon statistics is considered [30]. As illustrated in the reference, a uniform relation of statistical characteristics between 
the localized phonon and electron current does not exist. The relation is determined by the system parameters. As an example, in 
the parameter area where the charge-oscillator coupling is weak and two tunneling barriers are very asymmetric, one can expect 
a super-Poissonian current Fano factor associated with sub-Poissonian phononic population. The sub-Poissonian statistics of 
phonon distribution is also predicted in a system that a resonator coupled to a superconducting single-electron transistor J31I1 . 
They found the system behaves as the micromaser and can generate number-squeezed state of the resonator. Since the phonon 
number distribution would be narrowed in the squeezed state B23I1 . the phonon noise is reduced to be sub-Poissonian . 

Finally, we conclude that a master equation is developed to describe coherent dynamics of the electromechanical quantum 
transport as a function of the bias voltage and temperature. In the Born-Markovian approximation, the general master equation 
can be explicitly derived including position dependence of tunneling rates and Fermi distribution functions of the electron 
reservoirs. The mechanical motion of the charged QD is modeled by the quantum harmonic oscillator and coherently coupled 
to the electron transfer. The equation of motion is numerically solved in the Fock state Hilbert space of electron and phonon. 
It is found that the steps of current are related to the steps of mean energy in the mechanical oscillator. The dissipation rate 
of the resonator significantly affects the current intensity, especially when it is comparable with the bare tunneling rate. The 
gate voltage can be applied to control the current with discontinuous quantities. In the extremely low bias voltage, only the 
channel of the ground state is opened. In this case, we observe sub-Poissonian noise of electron transfer. When the bias is high 
enough to drive the oscillator partially into its excited states the system current manifests super-Poissonian noise. It reveals, in 
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the Coulomb blockade regime, positive correlation is generated between the multi-channel electron transports. If the dissipation 
rate of the vibrational modes is much faster than bare tunneling rates, the channels of excited states become unimportant and 
the shot noise would be suppressed remarkably. In higher temperature, smoothing of the current steps is observed. This event 
is almost accompanied by the increase of current shot noise due to temperature. Our present work enriches the researches on 
phonon assisted electron tunneling and the nano electromechanical oscillator. 
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